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ABSTRACT 

We evaluate the broadband afterglow emission of very high redshift gamma-ray bursts 
(GRBs) using standard relativistic blastwave models with both forward and reverse 
shock components. For a broad range of parameters, a generic property for GRBs at 
redshifts z ~ 5-30 is that the emission peaks in the millimeter to far-infrared bands 
with milli-Jansky flux levels, first at a few hours after the burst due to the reverse 
shock, and then again for several days afterwards with somewhat lower flux due to 
the forward shock. The radio, submillimeter and infrared continuum emission should 
be readily detectable out to z ^ 30 by the Atacama Large Millimeter Array (ALMA), 
Extended Very Large Array (EVLA), Square Kilometer Array (SKA) and other facili- 
ties. For relatively bright bursts, spectroscopic measurements of molecular and atomic 
absorption lines due to ambient protostellar gas may be possible. Utilizing models 
of primordial protostellar clouds, we show that under certain conditions, appreciable 
absorption may be caused by HD rotational transitions even in metal-free environ- 
ments. After sufficient metal enrichment, absorption from CO rotational transitions 
and [OI] fine-structure transitions can also become strong. With appropriate observ- 
ing strategies in combination with optical telescopes, ALMA and/or SKA may be able 
to detect such lines, offering a unique probe of physical conditions in individual Pop 
III and early Pop II star forming regions. We also remark on potential near-infrared 
absorption features due to electronic transitions of H2. 

Key words: radiation mechanisms: non-thermal - gamma-rays: bursts - ISM: 
molecules - galaxies: high-redshift - submillimctre - radio lines: general 



1 INTRODUCTION 

Gamma-ray bursts (GRBs) hold the promise of being pow- 
erful probes of the universe at very high redshifts* (Loeb 
2003, Djorgovski 2004). Their presumably close link with 
the collapse of massive stars indicates that they may occur 
up to the earliest epochs of cosmic star formation (Bromm 
& Loeb 2002, Schneider, Guetta & Ferrara 2002), and their 
extreme luminosities in broad wavebands should allow them 
to be readily detected out to the highest relevant redshifts 
(Lamb & Reichart 2000, Ciardi & Loeb 2000, Gou et al. 
2004), a realistic prospect with the launch of the SWIFT 
satellite mission (Gehrels et al. 2004). 

The formation era of the first stars in the universe is 
of considerable interest from many perspectives (Barkana & 

* In this paper, 'very high redshift' refers to z 10. 



Loeb 2001, Bromm & Larson 2004, Ciardi & Ferrara 2005, 
Glover 2005 and references therein). First, the pertinent 
physics in metal/dust-free conditions should be markedly 
different from the present day situation, most probably lead- 
ing to predominant formation of very massive, Population 
III stars. Second, such stars and their associated super- 
novae are expected to induce crucial feedback effects on their 
surroundings, including reionization and metal enrichment 
of the intergalactic medium, which will then strongly af- 
fect subsequent star and galaxy formation. After sufficient 
buildup of metals and dust, a transition to the more famil- 
iar mode of Population II star formation should take place. 
Finally, all these processes may be intimately related to the 
early formation and evolution of supermassive black holes 
and quasars. 

However, most of these inferences are theoretical and 
little is known with certainty, as observational clues are 
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very scarce. Aside from indirect indications such as evi- 
dence of early reionization in cosmic microwave background 
anisotropics (Spergel et al. 2003, 2006) or possible contri- 
butions to the cosmic near-infrared background (Santos, 
Bromm & Kamionkowski 2002, Salvaterra & Ferrara 2003; 
but see also Salvaterra & Ferrara 2006), detailed observa- 
tions of stellar populations in the corresponding redshift 
range z ~ 10 — 30 are still beyond the capability of cur- 
rent observational facilities. 

One of the most direct observational signatures of the 
star formation process is molecular line emission from the 
cooling, collapsing protostellar material, which generally 
arise in the radio, submillimeter to infrared wavebands and 
provide a powerful probe of star formation in the low redshift 
universe. In fact, emission lines from CO have been observed 
from galaxies with redshifts as high as z ~ 6.4 (Walter et 
al. 2003) , and next generation telescopes should certainly go 
deeper (Blain, Carilli & Darling 2004). However, it is uncer- 
tain whether this would be feasible much beyond z ~ 10, 
since 1) in the standard cosmological picture of hierarchical 
structure formation, the typical masses of collapsed, star- 
forming systems at the corresponding epochs are expected 
to be much smaller, and 2) the metallicity, and hence the CO 
abundance, should be much lower at earlier epochs of star 
formation. Emission line fluxes from H2 have also been esti- 
mated for very high-z populations, but they are typically so 
faint as to be undetectable by all but the most sensitive fu- 
ture observational facilities (e.g. Omukai & Kitayama 2003, 
Mizusawa, Nishi & Omukai 2005 and references therein; see 
however Ciardi & Ferrara 2001). 

A different approach to probe cool, star-forming gas 
at high redshifts is through absorption lines in the spec- 
tra of bright background sources such as quasars. Up to 
now, molecular absorption lines have been successfully de- 
tected in only a handful of systems up to z ~ 3.4, either 
through electronic transitions of H2 in the rest-frame UV 
(e.g. Ledoux et al. 2003, Reimers et al. 2003), or rotational 
transitions of molecules such as CO in the radio (e.g. Wik- 
lind & Combes 1999, Curran et al. 2004). More and higher z 
detections should undoubtedly follow with future telescopes 
(Curran et al. 2004, Blain et al. 2004), but the usefulness of 
quasars for such purposes may become limited at very high 
z: as with galaxies, sufficiently bright quasars (those lying on 
the tail of the luminosity function) may become much rarer 
in a standard, hierarchical cosmology, if the close relation 
between the masses of black holes and their host galaxies 
seen in the nearby universe also holds at early epochs. 

For absorption line studies at very high z, GRBs may 
offer an important alternative to quasars. Being the end re- 
sult of stellar collapse events, their typical luminosities are 
not expected to decrease strongly with redshift like quasars. 
They could even be typically more luminous at high z if 
black holes forming in GRB events are more massive than at 
lower z (Meszaros & Rees 2003). * GRBs also possess other 
advantages for this purpose, such as the simpler, nonther- 
mal intrinsic spectra devoid of strong emission lines, and the 
weaker perturbative effects on their environments (Barkana 
& Loeb 2003). The main disadvantage is their transience; in 

t However, whether Pop III stars can produce GRBs robustly is 
still an open question (Hcger et al. 2003). 



order to obtain adequate absorption line spectra, they must 
be observed while they are bright enough with suitable tele- 
scope settings, which requires fast response to GRB alerts 
with appropriate preparations. The implications of observ- 
ing absorption in high-z GRB spectra due to intergalactic 
hydrogen or metals have been discussed for rest-frame opti- 
cal/UV lines (e.g. Miralda-Escude 1998, Lamb & Reichart 
2000, Oh 2002, Furlanetto & Loeb 2003, Barkana & Loeb 
2003), and radio 21-cm lines (Furlanetto & Loeb 2002, Ioka 
& Meszaros 2005). 

After this paper was submitted, SWIFT discovered 
GRB050904 (Cusumano et al. 2006), which was successfully 
followed up at various wavelengths (e.g. Watson et al. 2006, 
Haislip et al. 2006, Boer et al. 2006, Frail et al. 2006). Spec- 
troscopic observations by the Subaru telescope have identi- 
fied its redshift as z — 6.29 (Kawai et al. 2006), the high- 
est of any GRB to date, and have already provided new 
constraints on cosmic reionization (Totani et al. 2006). The 
event marks an important first step in the use of GRBs as 
probes of the high-z universe, and showcases the potential 
of future studies at even higher z. 

This work focuses on the prospects for observing molec- 
ular and atomic absorption lines due to ambient protostel- 
lar gas in the radio to IR spectra of high-z GRBs, and the 
implications thereof for probing star formation at very early 
epochs. We first discuss in §2 the broadband spectra of high- 
z GRBs expected in standard afterglow models, including 
both the forward and reverse shock components. The de- 
tails of the spectral model are described in §A. The con- 
tinuum detectability in the radio to infrared bands is then 
assessed for various observational facilities such as ALMA 
and SKA. In §3, we present model calculations of primordial 
star-forming clouds incorporating detailed chemical reaction 
networks, from which the absorption line optical depths for 
different atomic and molecular transitions are evaluated. §4 
is devoted to a discussion of the detectability of such ab- 
sorption lines in GRB afterglow spectra by telescopes in- 
cluding ALMA and SKA, and the consequent implications. 
We conclude in §5. The following cosmological parameters 
are adopted, consistent with WMAP results (Spergel et al. 
2003): n m = 0.27, Q A = 0.73 and h = 0.71. 



2 BROADBAND AFTERGLOW EMISSION OF 
HIGH REDSHIFT GRBS 

2.1 Model spectra and light curves 

The broadband afterglow emission accompanying GRBs can 
be described quite robustly as nonthermal emission from 
electrons accelerated in relativistic shocks, which occur as 
the GRB-generating outflow is decelerated by the ambient 
medium (for reviews, see van Paradijs, Kouveliotou & Wijers 
2000, Piran 2004, Hurley, Sari & Djorgovski 2004, Zhang & 
Meszaros 2004, Meszaros 2006). The nonthermal electrons 
generated at the forward shock (hereafter FS) in the exter- 
nal medium give rise to synchrotron emission which typically 
spans the radio to X-ray range and decays relatively slowly, 
being observable up to days or weeks after the burst for 
GRBs at typical redshifts (z ~ 1). We describe the FS com- 
ponent as presented in the Appendix of Inoue (2004), which 
generally follows standard afterglow model prescriptions but 
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also includes a treatment of synchrotron self-absorption at 
high ambient densities. The values of the principal param- 
eters for the FS are fiducially taken to be consistent with 
those obtained from broadband modeling of well-observed 
afterglows (e.g. Panaitescu & Kumar 2001, Yost et al. 2003): 
isotropic-equivalent kinetic energy of blastwave E — 10 53 
erg, ambient medium density n = 1cm , half angle of jet 
collimation 9j — 0.1 rad, fractions of post-shock energy im- 
parted to relativistic electrons e e ,/ = 0.1 and magnetic fields 
£ s,/ = 0.01, respectively, and spectral index of electron dis- 
tribution pf = 2.2. 

Also important is the reverse shock (hereafter RS) that 
runs through the GRB ejecta, which has a much shorter 
lifetime but can involve amounts of energy dissipation com- 
parable to the FS. The associated synchrotron emission ini- 
tially peaks at lower frequencies and can be much brighter 
than the FS, typically dominating the optical and/or in- 
frared bands for z ~ 1 GRBs. Thereafter, it should rapidly 
fade and shift down in frequency toward the radio, falling be- 
neath the FS component a few days after the burst. This RS 
component can provide a consistent explanation of the op- 
tical flashes and radio flares observed in a number of GRBs 
(e.g. Piran 2004, Meszaros 2006 and references therein). 
Note that GRB050904 at z = 6.3 also exhibited a bright 
optical flash that can be interpreted as RS emission (Boer 
et al. 2006). 

Although the expected properties of the RS emission 
have been discussed extensively with regard to the opti- 
cal and/or radio bands, very few works have explicitly ad- 
dressed the details of the spectrum in between these bands 
and the behavior of the spectral peak at the self-absorption 
frequency (e.g. Nakar & Piran 2004, McMahon, Kumar & 
Piran 2006). Therefore, here we employ our own formalism 
for evaluating the RS component, as detailed in §A. As a 
unique feature, synchrotron self-absorption is accounted for 
utilizing the appropriate absorption coefficient instead of the 
commonly used blackbody limit approximation (e.g. Sari & 
Piran 1999a, b, Kobayashi & Zhang 2003), which is likely 
to be less reliable due to ambiguities in the time-dependent 
area of the emitting surface. This can be important for our 
purposes, as self-absorption generally determines the peak 
frequency and flux of the RS. Here we only treat the time 
evolution after the reverse shock has crossed the GRB ejecta, 
since the earlier evolution during shock crossing is not of im- 
mediate interest for the purposes of this work. In addition 
to E and n, further parameters important for the RS and 
their fiducial values are: initial bulk Lorentz factor of GRB 
ejecta Fo = 100, intrinsic GRB duration T = 100 s, and the 
RS counterparts of e e ,/, es,/ and pf, respectively e e<r — 0.1, 
ebj = 0.01 and p r = 2.2. We will assume throughout that 
e e ,r = e e ,/ and p r = pf and simply denote these as e e and 
p, as this choice is physically plausible and consistent with 
the current (limited) observations of optical flashes and ra- 
dio flares. In the case that we take e_g,r = es,/, we likewise 
simply write es. The time-dependent, broadband afterglow 
spectra for our fiducial parameter burst at z = 1 can be seen 
in Fig.Al. 

We thus evaluate the spectra and light curves of very 
high-z GRBs. Figs.l (a) and (b) show the broadband spectra 
of our fiducial GRB occurring at different redshifts from z = 
1 to 30, seen at fixed post-burst observer times t = 4 hours 
and 3 days, respectively. Figs. 2 (a)-(c) show the light curves 
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Figure 1. (a) Broadband spectra of GRB afterglows at different 
z as labelled at fixed post-burst observer time t = 4 hr, for our 
fiducial GRB. Overlayed are 5 a continuum sensitivities of various 
observational facilities, assuming integration times 50 % of t. (b) 
Same as a), but at t = 3 days. 



at selected observing frequencies for GRBs at z —15 and 
30. As mentioned above, only the post-RS crossing portion 
(t > t x ) is displayed for the light curves; the fluxes at earlier 
t should always be less than the first plotted point. 

Within t ~ 1 day, the two different components can be 
clearly distinguished in the spectra (see also Fig.Al). The 
FS contributes a broad, broken power-law spectrum extend- 
ing from radio to X-ray frequencies, with time-dependent 
break frequencies due to self-absorption, electron injection 
and cooling. On the other hand, the RS spectrum spans 
a narrower range in the millimeter to far-infrared, but can 
be significantly brighter in these bands, consisting of a self- 
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Figure 2. Light curves of GRB afterglows at z = 15 and 30 
at fixed observing frequencies (a) 10 13 Hz, (b) 300 GHz and (c) 
10 GHz. Only the post-RS crossing part (i > t x ) is shown; the 
flux at t < t x is always less than at t = t x . Overlayed are 5 
a continuum sensitivities of observational facilities as labelled, 
assuming integration times 50 % of t. 



absorption peak and a cooling cutoff. After t ~ few days, 
the RS has mostly faded away and the FS dominates the 
whole spectrum. 

Two striking features are evident for the RS component. 
First, the peak frequency (v ~ 200 GHz at t = 4 hr) stays 
almost independent of z despite the large cosmological red- 
shifting. Second, the flux at this frequency is also almost 
z-independent at z J; 5 and remains remarkably bright, 
at milli-Jansky levels even from z — 30. In contrast, the 
FS component is progressively dimmer with z at nearly all 
frequencies. These properties follow from cosmological time 
dilation and the fast evolving nature of the RS emission; a 
fixed observer time t corresponds to an earlier rest-frame 
time at higher z, when the RS component is much more 
prominent and possesses a higher rest-frame peak frequency. 
We can see this more specifically from the formulae in §A 
appropriate for our fiducial parameter set, which is the thin 
shell, slow-cooling case with the break frequencies ordered 
as v m>r < v a ,r < v c ,r- From Eqs.A22 and A26-A28, the RS 



peak frequency due to self-absorption depends on t and z as 
v a ocr^l + z)- 1 -"" with a a ~ -1.54 + 3.2/0 + 4), which 
is v a oc i" 103 (l + 2) ' 03 for p = 2.2. The dependence close to 
t^ 1 for fixed z means that the effects of time dilation almost 
exactly compensate for frequency redshift when compared 
at fixed t, resulting in a near-constant observed peak fre- 
quency. This behavior is quite insensitive to the value of p 
(e.g. v a oc f- 105 (l + z) 05 for p = 2.5). Similarly, from A30, 
the RS peak flux at v a is f(v a ) oc t a t (1 + z)- 1 ~ a J D~ 2 with 
czf ~ —0.97 + 1.6(p — l)/(p + 4), where Dl is the luminos- 
ity distance. For p = 2.2, this is f(v a ) tx t~ i :i (I + z) 2 - 3 D^ 2 . 
The steep decay with t at fixed z implies a large flux increase 
with z at fixed t, counteracting the dimming with distance. 

The greatly enhanced significance of the RS emission for 
very high- 2 GRBs has also been discussed recently by Gou et 
al. (2004). However, they did not account for self-absorption 
and ascertain the peak frequency of the RS, as their main 
interest was in the near-infrared to X-ray domain. 

Although to a lesser degree than the RS, similar trends 
are exhibited by the FS emission as well, as discussed by 
Lamb & Reichart (2000) and Ciardi & Loeb (2000). For our 
fiducial burst, the FS is slow-cooling after t ~ 15 minutes, 
in which case the FS peak frequency v v evolves with time as 
oc t~ 3/2 at fixed z (e.g. Sari, Piran & Narayan 1998). This 
implies v v oc (1 + z) 1 ^ 2 when compared at fixed t, differing 
only by a factor ~ 2 in the range z = 5 — 30 (e.g. v ~ 250-570 
GHz at t = 3 days; Fig.l (b)). The flux at fixed observing 
frequency and observer time is also of interest. For v > v p , 
f v falls with t so that time-dilation leads to a higher flux 
at higher z for fixed t; on the other hand, the spectrum in 
this region falls with v, so frequency redshift (K-correction) 
gives the opposite effect. For v < v p , the converse is true 
because f v rises with t and rises with v. In either case, the 
effects of time-dilation and frequency redshift act against 
each other to produce a flux which does not change greatly 
with z at fixed v and t. This expectation has indeed been 
observationally confirmed by Frail et al. (2006) in the radio 
band, at least up to z=6.3. Such effects should be even more 
dramatic for the RS due to its faster t-evolution. 

We have studied comprehensively how the fiducial re- 
sults above are affected when the main parameters are var- 
ied within a plausible range, and the fact that the after- 
glow spectrum manifests a peak in the millimeter to infrared 
range hours to days after the burst was found to be a quite 
robust property. As a demonstration of this statement, we 
show in Fig. 3 the results when n is assumed to depend on z 
as n = (1 + z) 3 cm -3 , which is one possible way that the in- 
terstellar medium density of high-z galaxies may evolve and 
a distinct alternative to our fiducial choice of constant n (see 
Inoue 2004 for more on this issue). The rapid increase of n 
with z means that with the same es, the post-shock mag- 
netic fields (in this case for both FS and RS) grow larger. 
This in turn implies that both synchrotron self-absorption 
and synchrotron cooling become stronger with z, leading to 
sharply peaked spectra by suppressing the lower and higher 
frequency portions, respectively. Intriguingly enough, when 



t Note that this trend continues only up to a certain z above 
which the observer time t < t x , when one starts seeing the pre- 
RS crossing part of the light curve and the flux decreases more 
strongly. 
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been assumed. 



observed at fixed t, the peak frequency is again almost inde- 
pendent of z and lies in the submillimeter band at t ~ few 
hours. 

We mention that our spectral model of §A is based on 
the approximation of taking either the ultrarelativistic or 
Newtonian limit descriptions for the RS dynamics. A more 
accurate evaluation may give quantitatively different results 
in the trans-relativistic regime (e.g. Nakar & Piran 2004). 
However, our intention here is to provide illustrative model 
examples of broadband afterglow spectra for their potential 
use in probing the high-z universe, rather than in compar- 
isons with real observational data. More detailed spectral 
modeling will thus be deferred to future work. 

2.2 Detectability of continuum emission 

Noting the very high flux levels even at the highest z, we as- 
sess the detectability of very high-z afterglows with current 
and future telescopes operating in the radio, submillimeter 
and infrared bands, concentrating on the fiducial case. In 
Figs.l - 3, the 5 a continuum sensitivities are overlayed for 
the current Very Large Array ( VL A) § , the future Extended 
Very Large Array (EVLA) Atacama Large Millimeter Ar- 
ray (ALMA) II , and Square Kilometer Array (SKA) ** facil- 
ities, as well as the current Spitzer telescope TT j all assuming 
that the integration times are 50% of t. 

Particularly interesting is the millimeter to far-infrared 
range where the RS component falls. Although observations 
here are generally more difficult than in neighboring bands 

§ http://www.vla.nrao.edu/ 

" http://www.aoc.nrao.edu/evla/ 

II //www. eso.org/projects/alma/; Wootten (2002) 

** http://www.skatelescope.org; Carilli & Rawlings 2004 

TT http://ssc.spitzer.caltech.edu/ 



due to strong atmospheric absorption and other factors, the 
prominent, mjy- level RS contribution allows the emission 
to be readily detectable out to z ~ 30 by ALMA with in- 
tegration times of a few hours. Even after the RS decays at 
t 1 day, the light curves in these bands generally brighten 
again at t ~ several days, due to the passage of the FS 
peak frequency (Figs. 1 (b) and 2 (b)). Although the flux 
of this second FS peak in the light curve is about an order 
of magnitude fainter than the first RS peak, the longer al- 
lowed integration times should make its detection by ALMA 
also feasible. Note that the limited sensitivity of currently- 
operating telescopes in these wavebands have so far enabled 
only the brightest, low-z bursts to be detected (Smith et al. 
2005a, 2005b, Kuno et al. 2004, Kohno et al. 2005). 

In the radio domain below few tens of GHz, the RS is 
never really important and the flux is much less at early 
times. However, the radio flux gradually brightens to 0.01- 
0.1 mjy as the FS peak frequency reaches these frequencies 
at t ij 10 days (Fig. 2). This should be well within reach of 
the generally superior sensitivities of radio facilities. Even 
the current VLA, which has played a crucial role in elu- 
cidating the properties of radio afterglows for GRBs ob- 
served to date (e.g. Frail 2003, Berger 2003), may detect 
very high-z bursts as well, albeit with long integration times 
of a few days. EVLA should be capable of doing this in sev- 
eral hours and SKA possibly in a few minutes. These points 
were also discussed recently by Ioka & Meszaros (2004). In 
the mid- to near-infrared, Spitzer should also have a chance 
with ~hour to ~day integration times, not to mention other 
ground based telescopes and the future James Webb Space 
Telescope (JWST) at near-infrared wavelengths (Ciardi & 
Loeb 2000, Gou et al. 2004). 

The event rate of GRBs at high-z is presently rather 
uncertain. Extrapolating the redshift distribution of bursts 
observed by SWIFT so far, Jakobsson et al. (2006) conclude 
that the fraction of bursts occurring at z > 5 should be 
7 %-40 %, which is in agreement with recent theoretical 
expectations (e.g. Natarajan et al. 2005, Bromm & Loeb 
2006). Bromm & Loeb (2006) also predict that the rate at 
z > 10 could be up to several per year. Although a detailed 
assessment of the rate is beyond the scope of this paper, 
it is conceivable that at least several detectable GRBs per 
year occur at z J> 5, and some even beyond z ~ 10. Their 
unambiguous identification, however, requires dedicated ob- 
servational efforts (see Sec. 4. 2). 

Successful detections of very high-z GRB afterglows in 
the above bands would have important implications. First, 
they may simply be the highest redshift sources of any kind 
to be detectable at these frequencies. The very first stars 
in the universe are expected to appear at z Jl 30 (Miralda- 
Escude 2003), and these may end up in GRBs, providing 
the "first observable light" in the radio to far-infrared wave- 
bands. Second, the time-dilated evolution of very high-z 
GRBs may prove to be valuable for studying the physics of 
the GRB outflow itself, by enabling the RS emission com- 
ponent to be observed more comprehensively than may be 
possible at low z. Detecting the fast decaying RS signature 
for GRBs at typical distances (z ~ 1) in the optical band 
has proved to be challenging, with only a handful of suc- 
cessful cases (e.g. Akerlof et al. 2000, Roming et al. 2005). 
While a number of explanations are possible for the paucity 
of detections (Meszaros 2006 and references therein), a full 
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understanding of the RS phenomenon may not be achieved 
without covering the millimeter to infrared range, where a 
significant portion of the RS spectrum might lie (Fig.Al and 
1). Observational facilities in these bands may possibly have 
a better chance of detecting the RS component for higher z 
GRBs, thanks to time dilation which allows more time for 
telescope pointing. 

However, the intriguing possibility which we highlight 
in the following is that some GRB afterglows may be bright 
enough in the radio to infrared bands for high resolution 
spectroscopic observations of molecular and atomic absorp- 
tion lines due to ambient cold gas, which could have far- 
reaching implications for exploring the earliest epochs of star 
formation in the universe. 



3 PRIMORDIAL STAR-FORMING CLOUDS 
AND ABSORPTION LINE OPTICAL 
DEPTHS 

We now turn to discuss how much absorption may be ex- 
pected by molecular and atomic fine-structure lines in star- 
forming environments at very high redshifts, within which 
GRBs may arise. Our estimates are based on Omukai (2000, 
2001) and Omukai et al. (2005), who constructed models of 
cooling, gravitationally collapsing clouds with varying levels 
of metallicity and external irradiation. Although the dynam- 
ics was treated with a one-zone approximation, the thermal 
evolution and non-equilibrium chemistry were calculated in 
detail by following a large number of reactions among the 
atomic and molecular species of H, He, C and O as well 
as dust grains. The calculations here follow Omukai et al. 
(2005), who has improved upon Omukai (2000, 2001) by in- 
cluding the chemistry of deuterium, in particular that of the 
crucial molecule HD, as well as the effects of the cosmic mi- 
crowave background (CMB). The potential influence of the 
radiation from the GRB itself is not included here, but is 
discussed in §4.1. 

A dynamically collapsing, spherical cloud without rota- 
tion or magnetic fields should possess a central core region of 
roughly constant density p and size given by the local Jeans 
length Xj = nc s /(Gp) 1 ^ 2 , where c 3 is the sound speed (e.g. 
Omukai & Nishi 1998). Focusing on this core region, the 
gas collapse is assumed to proceed on the free fall timescale 
tff = (3n/32Gp) 1 ^ 2 (note that the gravity due to dark mat- 
ter is important only in the earliest stages of the collapse). 
The initial number density and temperature are chosen to be 
n = 0.1cm -3 and T = 300 K in all cases; however, the calcu- 
lated subsequent evolution for a given composition quickly 
becomes independent of these values, as discussed in Omukai 
(2000). Defining the concentration y x — n x /n,H of a partic- 
ular species x as the ratio of its number density n x to that 
of hydrogen nuclei uh, the abundances of He and D nuclei 
are fixed to their primordial values j/ho = 9.72 x 10~ 2 and 
yr> = 4 x 10 -5 . The metallicity Z is parameterized by its 
fraction with respect to the local interstellar value Zq, as- 
suming that the relative abundances of heavy elements do 
not change. The model of Pollack et al. (1994) is adopted 
for the dust grain properties, and the dust-to-gas ratio is 
assumed to be proportional to Z. The redshift of the proto- 
stellar cloud is fiducially taken to be z — 15; this determines 



the CMB temperature, although its effect is minor for clouds 
of zero or very low metallicity. 

Another important parameter is the initial value of the 
H2 concentration j/H2,o before cloud collapse. For the very 
first generation of stars forming in H2 cooling-dominated 
mini-halos with virial temperatures T v < 10 4 K, the pre- 
collapse H2 concentration should be close to the primordial 
intergalactic value j/H2,o — 10~ 6 , since gas falling onto mini- 
halos is only weakly shocked. However, for later generation 
stars forming in more massive, atomic cooling-dominated 
halos with T v > 10 4 K, the initial H2 concentration can 



be as high as 2/H2,o 



10" 



Here the gas should first 



be strongly shocked and ionized at the virialization shock, 
which then lead to effective H2 formation mediated by abun- 
dant free electrons left over from non-equilibrium recombi- 
nation in the post-shock layer (Shapiro & Kang 1987, Oh 
& Haiman 2002). Similar considerations apply to gas that 
has been shocked and compressed by supernova blastwaves 
from a preceding generation of stars (Ferrara 1998, Bromm, 
Yoshida & Hernquist 2003) . We consider below the two cases 
of 2/H2.0 = 10~ 6 and 10~ 3 . 

The optical depth to absorption for each molecular or 
atomic line transition from level i to j is evaluated by 

where Aij is the Einstein coefficient for spontaneous tran- 
sitions, is the transition frequency, gt is the statisti- 
cal weight of level i, vd = {2ksT / PxTUh) 1 ^ 2 is the veloc- 
ity dispersion for a species with molecular weight p x , and 
l s h = min(i? c , Asth) is the shielding length, taken to be the 
smaller of either the core radius R c = Aj/2 or the Sobolev 
length As t h = vd /(dv/dr) = Svotff. This gives an esti- 
mate for the line optical depth through the core region of 
the cloud, the precise value depending on the exact geometry 
and specific line-of-sight. 



3.1 Zero metallicity clouds 

We first consider a completely metal-free (Z — 0) cloud. 
Fig. 4 shows the optical depth r for the two lowest-lying ro- 
tational transitions of H2 (J = — ► 2, v = 28.3 p,m and 
J — 1 — ► 3, v = 17.0 Atm) and the lowest-lying rotational 
transition of HD (J = — ► 1, v = 112 /im) as a function 
of the total number density n in a Z — cloud at redshift 
z = 15. Two different initial H2 concentrations have been 
assumed, yn2,o = 10~ 6 and 10~ 3 , as discussed above. There 
is a one-to-one relation between the time after collapse and 
the density (and hence the temperature) ; a given value of n 
corresponds to a particular time during the cloud evolution 
when the core has a size Aj ~ 1.1 x 10 18 (T2/n6) 1,/2 cm and 
is collapsing on a timescale i/y ~ 5.2 x 10 4 ng 1//2 yr, where 
the notations n = 10 6 n6cm~ 3 and T = WOT2K have been 
used. The evolutionary trajectories in the n — T plane for 
selected cases can be found in Omukai et al. (2005). 

We see that the optical depths for H2 in all situa- 
tions are generally quite low unless very high densities are 
reached. That for HD is similarly low if £/H2,o = 10~ 6 . Thus, 
for the very first generation protostellar clouds in mini-halos, 
strong absorption lines are not expected unless the line of 
sight passes through the small, innermost region of the col- 
lapsing core. 
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Figure 4. Absorption line optical depth r versus number density 
n for the rotational transitions J = — > 2 and J = 1 — * 3 of 
H2 and J = — ► 1 of HD, in protostellar clouds with 2 = 15, 
metallicity Z = and initial H2 concentrations J/H2.0 = 10~ 6 
and 10 -3 . 
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Figure 5. Absorption line optical depth r versus number density 
n for the rotational transitions J = — » 2 of H2 and J = — » 1 
of HD, in protostellar clouds with 2 = 15, J/H2,0 = 10 -3 and 
metallicities Z/Zq = 0, 10~ 4 and 10~ 2 . 



Contrastingly, in the case of yn2,o = 10~ 3 , the HD 
optical depth manifests an interesting behavior, increasing 
rapidly with n and approaching r ~ 1 at relatively low den- 
sities n ~ 10 6 — 10 7 cm -3 , and then declining temporarily 
at higher n. Initially, the temperature keeps decreasing and 
becomes low enough (T Sa 150K) for a large fraction of the 
available D to form HD, whereas the H molecular fraction 
still remains at ~ 10 -3 . HD also has a much larger A coeffi- 
cient thanks to its dipole moment, which H2 lacks. These two 
facts allow the HD optical depth to exceed that for H2 by 
a few orders of magnitude. However, HD reaches local ther- 
modynamic equilibrium above n ~ 10 5 cm -3 , and thereafter 
the temperature keeps rising until much of the HD becomes 
dissociated, leading to the drop in r. Note that around this 
temperature minimum, cooling is dominated by HD and al- 
lows the gas to reach lower temperatures (T ~ 50K) than 
is possible with H2 alone. 

Thus, for later generation star formation in relatively 
massive halos or supernova shells where j/H2 = 10 -3 can be 
realized, HD may cause sufficiently strong absorption even in 
metal-free environments, potentially offering a unique probe 
of Population III star formation (§4.3). 



3.2 Low metallicity clouds 

The line optical depths for clouds with ym.o = 10~ 3 and dif- 
ferent levels of metallicity are shown in Figs. 5 and 6. Higher 
Z allows more efficient H2 formation through dust grain sur- 
face reactions and consequent cooling to lower temperatures. 
This also induces HD formation, which cools the gas even 
further by becoming the dominant coolant for some interme- 
diate range of densities depending on Z (Omukai et al. 2005) . 
Unlike the Z — case, HD is not immediately destroyed at 
n <; 10 5 cm -3 because the temperature can remain low. All 




Figure 6. Absorption line optical depth r as a function of number 
density n for the rotational transitions J = — * I, J = A — >5 and 
J = 9 — > 10 of CO, and the lowest-lying fine structure transition 
of [OI], in protostellar clouds with 2 = 15, yu2 = 10 — 3 and 
Z/Z Q = 10~ 4 and 10" 2 . 



these facts contribute to a generally higher r with higher Z 
for both H2 and HD, as can be seen in Fig. 5. 

Fig. 6 displays the optical depths for different rotational 
transitions of the CO molecule (J = — > 1, v = 2.6 mm; 
J = 4 -> 5, v = 520 ^m and J = 9 -> 10, v = 260 fim). 
When Z = W~ 2 Zq, the optical depths of CO lines can ex- 
ceed unity at densities as low as n ~ 10 3 -10 5 cm . Even for 
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a very low metallicity of Z — 10~ 4 Zq, the higher-excitation 
transitions of CO can become strong at n ^ 10 7 cm -3 . 
Since the critical metallicity for transition from a Pop III, 
high mass dominated mode of star formation to a Pop II, 
low mass dominated one is expected to be in the range 
Z C rit/Z Q = 1CT 6 -I(r 4 (Bromm et al. 2001, Schneider et al. 
2002, 2003, Omukai et al. 2005), CO lines should constitute 
a valuable probe of the very early Pop II era, and possibly 
of the crucial transition epoch as well. We also see that for a 
given Z, the relative strengths of different rotational transi- 
tions vary with n (and hence T) in interesting ways. If more 
than one CO line can be measured for the same cloud, the 
line ratios should provide a useful diagnostic of the cloud 
physical properties (§4.3). 

Also shown in Fig. 6 is the optical depth to the lowest- 
lying atomic fine structure transition of [OI] (y = 63.2 lira), 
which approaches ~ 1 at n ~ 10 4 cm~ 3 if Z = \Q~ 2 Zq. 
That for [CII] is about an order of magnitude less. These 
absorption lines should be important for probing the earliest 
stages of the collapse, as [OI] and [CII] are thought to be 
the main cooling agents at this time (§4.3). 

We mention that the one-zone models presented here 
are intended only to provide order of magnitude estimates of 
the absorption line optical depths in primordial protostellar 
clouds. More realistic evaluations should be carried out on 
the basis of 3-D simulations (e.g. Abel, Bryan & Norman 
2002, Bromm, Coppi & Larson 2002, Yoshida et al. 2006, 
Gao et al. 2007), if the chemistry of HD or metals/dust can 
be fully incorporated within them in the future. 



4 OBSERVABILITY OF ABSORPTION LINES 
AND IMPLICATIONS 

In §2, we saw through standard relativistic blastwave mod- 
els that very high-2 GRB afterglows should have spectra 
robustly peaking around the submillimeter band, initially 
at a few hours after the burst with milli-Jansky flux levels 
due to the reverse shock, and then for a second time several 
days afterwards with somewhat less flux due to the forward 
shock. In §3, model estimates were presented for absorption 
line optical depths in primordial protostellar clouds which 
may host very high-z GRBs, demonstrating potentially ap- 
preciable absorption due to HD even in metal-free g 
well as strong absorption from CO and [OI] in the case of 
sufficiently metal-enriched gas. Bringing together these dis- 
cussions, here we assess the observability of molecular and 
atomic absorption lines due to cool, protostellar gas in very 
high-z GRB afterglows and its implications. 



4.1 Molecular gas in GRB environments 

The first important issue to be addressed is how much 
molecular gas may actually exist in the vicinity of GRBs 
when they occur. GRBs likely result from core collapse at 
the endpoint of a massive star's life, and these progenitor 
stars should be born embedded in the dense centers of star- 
forming clouds like those discussed in §3. However, a massive 
star emits copious amounts of UV radiation during its life- 
time, which is very effective in photodissociating molecules 
within the star's vicinity, before the GRB goes off. In fact, 



radiation from a single massive star may be enough to de- 
stroy the entire molecular content of a host mini-halo with 
T v < 10 4 K, at least from simple considerations assuming 
spherical symmetry (Omukai & Nishi 1999, Glover & Brand 
2001; see also Johnson, Greif & Bromm 2007 for recent 3-D 
simulations). Radiation from the GRB itself can also con- 
tribute to photodissociation; however, the total number of 
UV photons is at most ~ 10 61 even for a bright UV flash 
(e.g. Draine & Hao 2002), which is less than that expected 
from the progenitor star (e.g. Ciardi et al. 2001). Gamma- 
rays and X-rays are even less important because of the fewer 
number of photons and smaller absorption cross sections. 
At any rate, for mini-halos with T v < 10 4 K, the line opti- 
cal depths should be small anyway because of the low yn2,o 
(§3.1), so it may be unlikely that any absorption lines would 
be observable for GRBs emerging from the very first gener- 
ation of stars in the universe. 

The prospects may be better for GRBs from later gen- 
eration stars forming in halos with T v > 10 4 K. First, larger 
absorption line depths are realized due to the higher yH2,o 
(§3.1), as well as the larger gas masses in such halos. Second, 
the gas in these large halos can collapse into numerous clouds 
which then fragment into multiple clumps, each of which 
can become a star-forming core (e.g. Oh & Haiman 2002, 
Bromm, Coppi & Larson 2002). Such situations are reminis- 
cent of massive star formation in giant molecular clouds in 
the Galaxy today (e.g. Lada & Lada 2003). Within a clump 
hosting a massive star, destruction of molecules is likely to 
be effective due to photodissociation as well as heating by 
photoionization and/or stellar winds. However, other clumps 
within the same cloud that are sufficiently far away and have 
yet to form a star inside should be free from such adverse 
effects, because the molecular column densities are large 
enough for these clumps to be self-shielded from external 
UV radiation. For the example of our zero-metallicity cloud 
with 2/H2,o = 10~ 3 , regions that approach n ~ 10 6 cm~ 3 
and r ~ 1 for the HD line correspond to column densi- 
ties TVhd ~ 10 18 cm~ 2 , which far exceed the minimum col- 
umn density required for self-shielding, TVhd ~ 10 13 cm~ 2 
(Draine & Bertoldi 1996). Thus, such pristine clumps lying 
in the foreground of a GRB can potentially cause interest- 
ing absorption lines of HD. Note, however, that the radius 
~ 0.3 pc and lifetime ~ 5 x 10 4 yr expected for the high r 
core (§3.1) may be much smaller compared to those for the 
whole clump, estimated to be ~ 10 pc and 2 x 10 6 yr for an 
initial clump density n ~ 10 3 cm~ 3 as appropriate in halos 
with T v > 10 4 K (Oh & Haiman 2002). The probability of 
these clump cores occupying the GRB line of sight depends 
on their uncertain spatial and temporal distribution and is 
presently difficult to predict, but quantitative and realistic 
estimates should be tractable in the future through 3-D sim- 
ulations. 

A further possibility arises if the photodissociation re- 
gion (PDR) around a massive star can sweep up the sur- 
rounding material and trigger further star formation. Such 
situations have indeed been observed inside our Galaxy (De- 
harveng et al. 2003, Hosokawa & Inutsuka 2005), and may 
also be possible in very high-z environments, although this 
is yet to be fully investigated. Newly collapsing clumps on 
the PDR periphery would have suffered little photodissocia- 
tion, and a relatively high j/H2,o is plausible due to the shock 
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compression. The chances may be better for such clumps to 
induce HD absorption lines. 

After metal enrichment has proceeded to above a crit- 
ical value Zcrit/ZQ = 1CF 6 — 1CT 4 , cooling by metals and 
dust becomes important. Besides the higher collapse densi- 
ties attained through more efficient cooling, the buildup of 
dust also leads to more effective formation of H2 on dust 
grain surfaces, as well as greater shielding against external 
UV radiation. For example, in our case with Z/Zq = 1CP 2 , 
regions with n J; 10 cm~ and r <; 1 for the lower-excitation 
CO lines have H2 column densities Nh2 ~ 10 21 cm -2 , which 
is sufficient for shielding of CO, not to mention H2 and HD 
(Lee et al. 1996). The expected optical depths can be high 
through most of the collapsing cloud, especially for CO lines, 
those for HD are also higher compared to when Z = (§3.2). 
It is therefore quite probable that strong molecular absorp- 
tion lines are formed in such metal-enriched environments. 

While having no destructive capability, the radio to in- 
frared emission from the GRB can potentially alter the level 
populations of ambient molecules through radiative excita- 
tion. Its importance can be approximately gauged through 
the simple example of a two-level system, described in de- 
tail in Appendix B of Stahler & Palla (2004). We assume 
a homogeneous medium of density n composed of particles 
with two levels separated by energy hi>o, immersed in an 
isotropic radiation field of intensity I in steady state con- 
ditions. Radiative excitation plays a significant role when 
n n' cr = n cr (1 + c? I /21%Vq), where n cr = Ao/70 is the crit- 
ical density above which collisions dominate and the system 
attains LTE in the absence of an external radiation field, be- 
ing expressible in terms of the Einstein coefficient for spon- 
taneous emission Ao and the coefficient for collisional deex- 
citation 70 between the two levels. In all cases of our interest, 
n ^> n cr . We may define a critical radiation intensity above 
which radiative excitation is important, 

_ 2hv% n 

J-cr — 2 ' \ Z ) 

to be compared with the mean intensity Iqrb of our fiducial 
afterglow at a distance d from the GRB in the molecular 
cloud. 

For the conditions relevant to the J = — > 
1 HD line (f =2.68 THz) in our Z = 0, jy H 2,o = 
10~ 3 case (§3.1), n/n cr ~ 30 and I cr ~ 8.5 x 
10 _9 erg s _1 cm _2 Hz _1 sr _1 . When the afterglow is bright- 
est at this rest-frame frequency at t ~ 10 min (Fig.Al), 
Iqrb ~ 10 _9 erg s _1 cm _2 Hz _1 sr _1 (ci/10pc) _2 , which drops 
thereafter by 1 to 2 orders of magnitude hours to days af- 
ter the burst. Therefore, the afterglow radiation will have 
no bearing on the excitation states of HD unless they lie 
at d <C 10 pc. As discussed above, the vicinity of the GRB 
should be devoid of molecules anyway due to the actions 
of the progenitor star. The effect could be more relevant 
for some of the CO lines in metal-enriched clouds (§3.2). 
For example, where t J5 1 for the J = — > 1 CO line 
(^o=H5 GHz) in the Z/Zq = 10~ 2 case, n/n cr ~ 100 
and I cr ~ 2.2 x 10~ 12 erg s _1 cm~ 2 Hz _1 sr _1 . At the cor- 
responding frequency, the fiducial afterglow has Iqrb ~ 
10 _10 erg s _1 cm _2 Hz _1 sr _1 (c?/10pc) _2 hours to days after 
the burst, so CO molecules within d ~ 100 pc may have 
their lower excitation states affected. A realistic evaluation 
of such effects is outside the scope of this work and set aside 



for future studies, requiring spatially-dependent calculations 
with time-varying, anisotropic radiation fields. We proceed 
below under the assumption that the molecules are located 
at d 10 — 100 pc where radiation excitation is negligible, in 
line with the above discussion of multiple foreground clumps 
in large clouds and/or clumps precipitated by PDR shocks. 

4.2 Detectability of absorption lines 

Given that the circumstances envisioned above are realized 
for at least some GRBs, the detectability of the absorption 
lines with different observational facilities can be appraised. 
The absorption line widths should be dominated by velocity 
gradients from the bulk infall motion in the collapsing cloud, 
which is roughly 

Av ~ Xj/tff ~ 6.7 km/s T 2 1/2 . (3) 

Since the temperatures are generally T ~ 30 — 100 K when 
large r values occur, a velocity resolution of at least 3 km/s 
would be necessary to detect such lines. We note that dif- 
fuse regions in which turbulent motions are important are 
unlikely to have sufficient optical depth for absorption lines. 

In Fig. 7, we show the broadband spectra of a relatively 
bright burst with E — 10 54 erg (all other parameters being 
fiducial), at different redshifts from z —1 to 30 and post- 
burst observer times t = 1 and 10 days. Marked on each 
spectra are the redshifted frequencies for the lowest-lying 
transitions of H2, HD, CO and [OI]. To be compared are the 
5-cr sensitivities of ALMA, EVLA and SKA, for both con- 
tinuum and spectroscopic observations at 3 km/s resolution, 
assuming integration times 50 % of t. Note that this spec- 
troscopic sensitivity is for optically thick absorption lines; 
to detect lines with small or moderate optical depths at the 
same resolution and significance, the required sensitivity (or 
alternatively, the required background continuum flux) is 
higher by a factor (1 — e~ r ) _1 (~ r _1 for r <C 1). Fig.8 dis- 
plays the light curves of GRBs at z =15 with different values 
of E and n, at observing frequencies u=10 GHz (near the 
low-excitation CO transitions) and 230 GHz (near the HD 
or [OI] transitions), in comparison with the spectroscopic 
sensitivities of the aforementioned telescopes. 

The redshifted absorption lines for the lower-excitation 
transitions of CO are mostly in the radio regime, while those 
for HD and [OI] , as well as the higher-excitation transitions 
of CO all fall in the submillimeter band. There are inter- 
esting prospects for all of these lines to be detectable. For 
GRBs with sufficiently large E J> 4 x 10 53 erg, ALMA may in 
fact have two opportunities to detect the submillimeter lines, 
first due to the passage of the bright RS peak at t 1 day, 
and second due to the passage of the FS peak around t ~ 
10 days (however, see below on practical difficulties with the 
first RS peak). Note that large isotropic equivalent energies 
E <, 10 54 erg are not uncommon in observed GRBs at high 
redshifts (Amati 2004), including GRB050904 at z = 6.3 
(Cusumano et al. 2006). " A burst with average E but 
relatively high n ii 10 2 cm~ 3 should also manifest a bright 
enough FS component at t ~ 1 day (Fig.8). 

++ The true total energy may be less due to collimation, but 
collimation effects in the afterglow should not be important in 
very high- z afterglows until t <; 10 days. 
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Figure 7. Broadband spectra of GRB afterglows with fiducial 
parameters except for E = 10 54 erg, at different z as labelled 
and fixed post-burst observer time (a) t = 1 day and (b) 10 days. 
The redshifted frequencies for the lowest-lying transitions of H2 
(crosses), HD (circles), CO (triangles) and [OI] (vertical bars), as 
well as the Ly a (diamonds) and HI 21cm (squares) transitions 
are indicated on each spectra. Overlayed are 5 cr continuum and 
3 km/s resolution spectroscopic sensitivities of various observa- 
tional facilities, assuming integration times 50 % of t. 



In the radio, the chances are even better. The CO lines 
should be detectable by EVLA for GRBs with £^4x 10 53 
erg at t ^ 10 days. SKA, with its outstanding sensitivity, 
should do the same for typical energy bursts in less than a 
day. With longer integration times, SKA may even be capa- 
ble of measuring weak lines with r 0.1, as well as strong 
lines in subluminous bursts. 

Besides adequate sensitivity and spectral resolution, 
some important technical constraints pertaining to ra- 
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Figure 8. Light curves of GRB afterglows with different E and 
n (as labelled by their logarithmic values), at 2 = 15 and at fixed 
observing frequencies (a) 230 GHz and (b) 10 GHz (only the post- 
RS crossing part t > t x is shown). Overlayed are 5 <r, 3 km/s 
resolution spectroscopic sensitivities of observational facilities as 
labelled, assuming integration times 50 % of t. 



dio/submm observations must be spelled out here. One is 
that interferometers such as ALMA or SKA require suffi- 
ciently bright calibrator sources near the intended targets 
in the sky in order to carry out absolute flux measurements. 
Depending on the actual distribution of suitable calibrators, 
which will only be known after ALMA begins operation, 
some restrictions may be imposed on regions of the sky in 
which GRBs can be observed. 

However, the most important caveat regards frequency 
coverage. For radio and submillimeter observations, achiev- 
ing high frequency resolution as is necessary here mandates 
that the instantaneous bandwidth be made correspondingly 
narrow. For example for ALMA, the number of spectral 
channels that the correlator can provide for an analyzing 
bandwidth B is 128(.B/2GHz) _1 per baseline for two po- 
larizations. A velocity resolution of 3 km/s at observing 
frequency v=230 GHz implies that the maximum allowed 
bandwidth would be 4B ~ 3.2 GHz (T. Wilson, private 
communication), so at a given time, only a small portion 
of the spectrum can be covered. This is not a great con- 
cern when searching for lines in steady sources, as a broad 
spectral range can be scanned over time. However, transient 
sources such as GRBs pose a problem: even if one aims for a 
particular line with a definite rest-frame frequency, its red- 
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shift must be known beforehand to good accuracy (for the 
above example, to within 1 %) in order to begin observing 
at the appropriate frequency while the source is still bright. 

Thus a crucial prerequisite is that good redshift infor- 
mation for the GRB is quickly available by other means. 
Spectrographs on moderate-sized optical telescopes are cur- 
rently indeed able to measure redshifts for GRB afterglows 
through metal absorption lines with accuracies better than 
fractions of a percent within a day or so after the burst, 
at least up to z ~ 6 (e.g. Fynbo et al. 2006 and references 
therein). Doing so for bursts at higher z is complicated by 
contamination from infrared airglow, but can be achieved 
through spectroscopy of the Ly break by large telescopes 
with near-infrared capabilities. For example, spectrographs 
similar to the CISCO instrument on the Subaru telescope 
can obtain a J-band spectrum (covering the Ly break at 
z = 10) at resolving power R=200 with S/N=10 in 1 hour 
for a 3 /iJy source (T. Yamada, private communication), 
which is more than sufficient for the afterglow fluxes dis- 
cussed here (Fig. 7). 

Rapid dissemination of such redshift information over 
the GRB Coordinate Network is already being carried out 
extensively. In the ALMA era, successor missions to SWIFT 
such as EXIST or ECLAIRs ^ should be in operation. 
In synergy with gamma-ray and optical/infrared telescopes, 
successful detections of radio/submm absorpton lines in very 
high-z GRB afterglows should be feasible. For this purpose, 
the FS shock peak at t <; several days should be much more 
suitable than the early decaying RS shock peak, because 1) 
more time is allowed before the redshift is measured and 
the appropriate radio/submm observations can begin, and 
2) the FS emission should be more robust compared to the 
RS emission, which can be weak for some bursts with strong 
RS magnetic fields. 

In addition to the above lines, we have also indicated 
the 21 cm hyper-fine transition frequency (1.4 GHz) in the 
spectra of Fig. 7. Although strong 21 cm absorption may po- 
tentially be caused by host galactic disks at high-z (Furlan- 
etto & Loeb 2003), it can be readily seen that even with 
the sensitivity of SKA and a bright GRB, detecting such 
absorption lines in GRB afterglows should be very difficult 
unless z <, 10. This is because the redshifted 21 cm line lies 
deep in the low frequency, self-absorbed part of the afterglow 
spectrum where the flux is at micro- Jansky levels. These in- 
ferences were mentioned in Inoue (2004) and also agree with 
Ioka & Meszaros (2004). For studying 21 cm absorption at 
very high-z, radio-loud quasars may thus be more promising 
sources (Carilli, Gnedin & Owen 2002, Haiman, Quataert & 
Bower 2004), although the existence and number of such 
objects at z J> 10 is quite uncertain at this moment. 

4.3 Implications of detecting absorption lines 

For each kind of absorption line discussed above, there are 
tremendous implications if they can actually be detected. 

As mentioned in §3.1, HD can play a crucial role in the 
cooling of zero or low metallicity gas in addition to H2. In 
some cases it can actually become the dominant coolant, 

§§ http://exist.gsfc.nasa.gov/ 

http:/ /www. oamp.fr/ECLAIRS/ 



cooling the gas down to few 10s of K and possibly even 
allowing the formation of relatively low mass stars (Uehara 
& Inutsuka 2000, Nakamura & Umemura 2002, Flower 2002, 
Nagakura & Omukai 2005, Johnson & Bromm 2006). The 
true importance of HD cooling in primordial star formation 
is not yet clear , but observations of HD absorption lines in 
high-z GRBs by ALMA should provide a quantitative test. 
If it can be demonstrated at the same time that the metal 
content of the gas is so low as to be insignificant for cooling 
(either through metal lines in the near- infrared or CO lines), 
it will offer a very unique signature of Population III star 
formation (albeit not for the very first generation of stars 
in mini-halos, but for later generation stars in more massive 
halos or triggered by supernovae/PDR shocks). 

It was shown that some high-excitation CO lines (e.g. 
J — A — > 5 to J = 9 ^ 10) can be relatively strong even 
when Z <; 10~ Zq (§3.2), which may be near the criti- 
cal metallicity for transition in star formation from high 
mass dominated, Population III to low mass dominated, 
Population II (Bromm et al. 2001, Schneider et al. 2002, 
2003, Omukai et al. 2005). This points to such lines (po- 
tentially observable by ALMA) being an important probe 
of the earliest phases of Pop II, and perhaps even of the 
critical transition era, which is of great interest from many 
viewpoints (e.g. Mackey, Bromm & Hernquist 2003, Scanna- 
pieco, Schneider & Ferrara 2003, Salvaterra & Ferrara 2003). 

At higher metallicity levels, the low-excitation CO lines 
(e.g. J = — ♦ 1 to J = 4 -> 5) can be very strong, be- 
ing readily observable by SKA and in some cases by EVLA. 
In fact, SKA should be sensitive enough to detect several 
different rotational transitions of CO from the same cloud 
during the time that the GRB afterglow is bright (Fig. 8 (b)). 
Measurement of line ratios will give information on physical 
conditions such as temperature and density for individual 
clouds, which may be impossible to do in any other way. 
Note that detecting the HD lines may also be possible for 
metal-enriched clouds and provide important additional in- 
formation. 

During the earliest stages of collapse, [OI] and [CII] are 
the main cooling agents for metal-enriched clouds (Omukai 
2000, Bromm & Loeb 2003, Omukai et al. 2005). Observing 
[OI] in absorption by ALMA will therefore constitute an 
interesting diagnostic of this phase, and may also offer a 
different perspective from CO lines on the Pop III to Pop II 
transition epoch. 



4.4 Near-infrared absorption by H2 

We briefly touch upon the possibility of observing electronic 
absorption features by H2 in the near-infrared. Absorption 
due to electronic transitions in the rest-frame UV are poten- 
tially very interesting because the A coefficients are much 
larger and the required column densities much smaller com- 
pared to ro-vibrational transitions. Indeed, as mentioned in 
§1, electronic absorption features have already been used 
for damped Lyman alpha systems at z <, 3.4 as an im- 
portant probe of H2, sensitive down to column densities 
Ah2 ^ 10 14 cm~ 2 (Ledoux et al. 2003, Reimers et al. 2003), 
and in one case even of HD (Varshalovich et al. 2001). How- 
ever, absorption from cold H2 in the Lyman- Werner bands 
(912-1110 A) will not be useful above z ~ 6, becoming com- 
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pletely obscured by Gunn-Peterson troughs shortward of Ly 
a from intergalactic HI (White et al. 2003). 

On the other hand, Hb in the local environment of GRBs 
may be more interesting. We have mentioned in Sec. 4.1 
that the bulk of the molecular gas in sufficiently distant 
clumps should be self-shielded from photodissociation by 
external UV radiation. However, the layers closest to the 
GRB can be affected by UV photons from the burst. Drainc 
(2000) and Draine & Hao (2002) have shown that if H 2 with 
Nh2 JJ. 10 18 cm -2 lies at distances R 10 19 cm around a 
GRB, which can indeed be the case for our clouds discussed 
above, the optical-UV emission from the GRB itself can 
pump the ambient H2 into vibrationally excited states. This 
will in turn induce strong, time-dependent absorption in the 
GRB spectrum up to rest-frame 1705 A, with characteristic 
features that can be identified through spectroscopy at reso- 
lution R J> 300. This may then provide an observable signa- 
ture of H2 even at z ^ 6 appearing in the near- infrared long- 
ward of Ly q, which we have also marked in Fig. 7. In fact, 
such a signature may have already been seen in GRB050904 
at z — 6.3 (Haislip et al. 2006). Although a detailed dis- 
cussion is beyond the scope of this paper, such electronic 
absorption features may offer interesting prospects for tele- 
scopes with fast response and near-infrared spectroscopic 
capabilities as a unique and valuable probe of H2 at very 
high-z. 



5 CONCLUSIONS 

The salient points of this work are summarized. We have 
evaluated the time-dependent, broadband afterglow spectra 
of very high-z GRBs using standard relativistic blastwave 
models with both forward and reverse shock components, 
paying particular attention to the latter. It was found that 
for a plausible range of physical parameters and external 
medium densities, GRBs in the redshift range z ~ 5 — 30 
exhibit spectra which robustly peak in the millimeter to in- 
frared range with milli-Jansky flux levels when observed a 
few hours after the burst, mainly due to the reverse shock. 
Following the decay of the reverse shock component, the 
light curve in these bands can brighten again for several 
days afterwards due to the forward shock. The continuum 
emission should be readily observable up to z ~ 30 and be- 
yond by telescopes such as ALMA, EVLA, SKA and other 
facilities in the radio, submillimeter and infrared. 

For relatively bright GRBs, high resolution spectro- 
scopic measurements of atomic and molecular absorption 
lines due to ambient protostellar gas may be possible. 
We have utilized models of zero- and low-metallicity star- 
forming clouds to show that under certain conditions, suf- 
ficiently strong absorption may be caused by transitions of 
HD, CO and [OI]. Such atoms and molecules should be im- 
mune to dissociation or excitation by radiation from the 
progenitor or the GRB as long as they lie 10-100 pc away 
from the GRB site. HD lines may constitute a unique probe 
of metal-free, Pop III star formation, albeit not for the very 
first generation of stars in mini-halos, but for later genera- 
tion stars forming in more massive halos or on the periph- 
eries of photodissociation regions. CO and [OI] lines can 
offer valuable information on the early Pop II era, possibly 
including the Pop III to Pop II transition epoch, and also 



provide diagnostics of physical conditions in the protostel- 
lar gas through line ratios. A cooperative observing strategy 
with gamma-ray and optical telescopes should enable the de- 
tection of HD, [OI] and high-excitation CO lines by ALMA, 
and low-excitation CO lines by SKA (and in some cases, 
EVLA) . Absorption features due to H2 in the near- infrared 
may also be a potentially interesting signature. 

Because of their transient nature, the detection of ra- 
dio/submm absorption lines in GRB afterglows may be chal- 
lenging in practice. However, they offer us a unique view into 
star- forming regions at very high redshifts (z J> 10), which 
may be quite difficult through emission lines or absorption 
lines in quasars for the reasons discussed in §1. On the other 
hand, at lower redshifts observationally accessible by emis- 
sion lines or quasar absorption lines, their steady fluxes give 
them an advantage. Lines in absorption and emission also 
provide different kinds of information, the former on local 
conditions along specific lines of sight, and the latter on 
globally integrated properties such as total mass. In these 
respects, different approaches using emission lines, absorp- 
tion lines in quasars and in GRBs should all be valuable in 
complementary ways for investigating early star formation. 
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APPENDIX A: BROADBAND AFTERGLOW 
EMISSION 

The forward shock (FS) component of the GRB afterglow 
emission is modeled as in the Appendix of Inoue (2004). This 
generally follows standard prescriptions for the synchrotron 
emission from electrons accelerated by a relativistic shock 
decelerating self-similarly in the ambient medium (e.g. Sari, 
Piran & Narayan 1998, Panaitescu & Kumar 2000), and in- 
cludes a treatment of synchrotron self-absorption at high 
ambient densities. As explained in §2, the main parameters 
are E, n, 9j, e e ,f, £bj and pf. The spectrum consists of 
power-law segments adjoining break frequencies at v m j and 
v c j, corresponding respectively to the injection energy ^m.j 
and cooling energy in the electron distribution, as well 
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as v a j, the self-absorption frequency. The peak flux f p j oc- 
curs at f p j = mm(v m j, v c ,f) when v a j < v P j, and at v a ,f 
reduced by self-absorption when v a j > v p j. The break fre- 
quencies and the peak flux evolve with time in characteristic 
ways through the evolution of the shock radius, bulk Lorentz 
factor and post-shock magnetic field. For more details and 
specific expressions, the reader is referred to Inoue (2004). 

For the reverse shock (RS), additional parameters are 
To, T, e e ,r, £B,r and p r , of which we take e e>r = e e ,/ = e e 
and p r = pf = p (§2). The dynamical behavior of the RS 
can be approximately classified into two cases, depending on 
the relation between the observed GRB duration T(l + z) 
which characterizes the initial ejecta width A = cT, and 
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the observed timescale for the forward shock to sweep up 1/T 
of the ejecta mass from the ambient medium (Sari & Piran 
1995, Kobayashi 2000). In the numerical expression, we have 
used the notation Q — 10 P Q P to represent physical quantity 
Q. When T(l + z) > tr, referred to as the thick shell case, 
the reverse shock becomes relativistic in the frame of the 
unshocked ejecta before the RS crosses the whole ejecta at 
t x ~ T(l + z), and the ejecta is decelerated to bulk Lorentz 
factor r x — F c , where 
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When T(l + z) < tr, referred to as the thin shell case, the 
reverse shock becomes marginally relativistic in the frame 
of the unshocked ejecta just when the RS finishes crossing 
the ejecta at t x — tr, and the ejecta is not significantly 
decelerated, T x ^ To- 

As with the FS, the RS spectrum is a broken power-law 
characterized by its injection frequency v m ,r and cooling fre- 
quency v c>r corresponding to break energies 7 m _ r and 7 c>r in 
the electron distribution, along with the self-absorption fre- 
quency v a ,r and peak flux f P;r . Its shape and evolution can 
be modeled in an approximate way by relating v m ,r, v c ,r 
and /p, r with their FS counterparts at the shock crossing 
time t — t x (Kobayashi & Zhang 2003, Zhang, Kobayashi 
& Meszaros 2003). From equality of the bulk Lorentz fac- 
tor and energy density in the post-shock regions across the 
contact discontinuity, we can derive the relations 
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where TZb = CB, r /eB,f, and Yf and Y r are the Compton pa- 
rameters for the FS and RS, respectively. As was assumed 
for the FS in Inoue (2004), Y r is approximated by its fast 
cooling value Y r = (y/ 4e e ,r/es,r + 1 — l)/2 at all times. Ad- 
ditionally, from energy conservation between the kinetic en- 
ergies in the initial ejecta and the shocked ambient medium, 



one gets the relation 
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Hereafter, all quantities pertain to the RS, and we drop 
the subscript r for brevity. After shock crossing (t > t y ), 
the time evolution of the above quantities can be obtained 
through the evolution of the shock radius R, bulk Lorentz 
factor T, number density Af, energy density £ and magnetic 
field B of the shocked ejecta as given in Kobayashi (2000). 
In the thick shell case, 
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Since 7 m oc £/M, v m oc r£>7^, and v c evolves in the same 
adiabatic way as v m after shock crossing (Sari & Piran 
1999a,b, Kobayashi 2000), 
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Similarly, as f p oc TBAfB 3 , 
f p oc t- 47/48 ~ t- (h98 . 
In the thin shell case, 
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Although the evolution before shock crossing (t < t x ) can 
also be deduced from suitable expressions for the evolution 
of T, etc. (e.g. Ioka & Meszaros 2004), it is outside of our 
interest here. 

The RS self-absorption frequency has often been esti- 
mated from the blackbody limit approximation (e.g. Sari & 
Piran 1999a,b, Kobayashi & Zhang 2003), which, however, 
suffers from uncertainties in the area of the emitting sur- 
face and its time-dependence. As a unique feature of this 
work, we utilize a more proper approach through the ap- 
propriate absorption coefficient. Denoting 7 P = min(7 m ,7 c ) 
and v p = min(^ m , v c ), we can use the self-absorption coeffi- 
cient averaged over an isotropic distribution of pitch angles 
(Granot, Piran & Sari 1999, Panaitescu & Kumar 2000) to 
express the self-absorption optical depth at frequency v p as 
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where f(p) = (p + 2)(p - l)/(3p + 2), V = 2 8 / 3 7r 5 / 2 /5r(5/6), 
and L is the comoving frame width of the shocked ejecta. 
The expressions for 7 m and 7 C are 

p — 2rn p To 



7m = e< 



and 



7c = 



p — 1 m e r 



l + z 



(A15) 



(A16) 



lQarc m e €BnT 3 t ' 

At t = t x , the following relations for Af and L hold approx- 
imately in either the thick and thin cases, 



My 



8Tj 
To 



(A17) 
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and 

L x -- 



r x ct x 

2(1 + 2) 



(A18) 



(Sari & Piran 1995, Kobayashi 2000). Thus, we can evaluate 
the self-absorption optical depth r PiX at v p at the shock 
crossing time by substituting t — t x and V — F x in Eqs.A14- 
A16 along with Eqs.A17 and A18. In numerical terms: 
when 7 C < j m (fast cooling), 

t p , x = 2.5 x 10- 7 /W(l + Yf^%E%\^Y-^T-n^) 
for the thick shell case, and 

r p , x = 1.5 x 10- 7 /(p)(1 + Y) 5 e 9 J 2 _ 2 Ei 3 n 7/2 r ,2 (A20) 

for the thin shell case; 
when 7 m < j c (slow cooling) , 



r p , x =0.34/(p) 



P - 



for the thick shell case, and 



o.9i/(p) 



e- 5 _ ie -^ 2 2 S 63 n- 1/2 r jT 2 -tA21) 



V2 t,l/3 l/6p-2/3 



-l c S,-2 iJ 53 



re 



■ 0,2 



(A22) 



for the thin shell case. 

The time evolution of t p after shock crossing can be 
derived in a similar way to Eqs.A8-A13, using B oc £ 1//2 , 
L cc Af~ 1 r~ 2 and r p oc J^LB^ 1 j p ''' . In the thick shell case, 



AT oc r 13 / 16 , L oc t 9/16 ,B oc t- 13/24 , 7P oc i^ 48 , 
so that 



whereas in the thin shell case, 



AToct- 6/7 ,Loct 16/3: 
so that 

.8/5 .1.6 
T v 0(.t — t 



,Bocr 4/7 ,7 P ocf~ 2/7 , 



(A23) 
(A24) 
(A25) 
(A26) 



Thus we have the value of t p at any t > t x . 

Writing v b — max(i/ m ,i/ c ), the self-absorption optical 
depth at an arbitrary frequency and time after shock cross- 
ing is then 

-5/3 

V < V p 

-(9+4) /2 

u p < v 4 A27) 

-(9+4)/2 

Vb < v 



T v = T p < 



-(p+5)/2 



I/f> 



The frequency at which t v = 1 gives the synchrotron self 
absorption frequency i/ a , 



i/p < 



3/5 
•p 

2/(9+4) 



2/(p+5) / ^b 



l-(q+4)/(p+5) 



Tp < 1 

r b < 1 < 

n > 1 



&28) 



where T b = T p (v b / 'i/ p ) _ ^ +4 - )/ ' 2 is the optical depth at v b . The 
case of T(, > 1 can be important for high ambient densities. 
The power-law index of the emission spectrum between 



v v and Vb is (q — l)/2, where q is the corresponding index 
in the electron distribution, being q — 2 for fast cooling and 
q = p for slow cooling. For the broadband spectrum /„ , three 
cases can be distinguished depending on the location of v a 
with respect to v v and v b : for v a < v p , 



1/3 



1/3 



"(9-l)/2 



V < V a 



Va < V < Vp 



(A2,9) 



-P/2 



"(9-l)/2 



Vp < V < Vb 



Vb < V 



for I/p < I/ a < I/;,, 



1/ 

l/ D / VI/, 
5/2 , 



i/p N 5/2 / 



"(9-l)/2 



(a 



"(9-l)/2 



(-) 



-p/2 



"(9-1)72 



"(9-l)/2 



V < V p 

V p < V < V a 

(A30) 

I/ n < V < V b 
Vb < V 



and for v b < v a , 



fp < 



v b 



5/2 



V a ) Wb 
(9"l)/2 



-p/2 



1/ < I/p 



(0 © 



-P/2 



?//, 



"(9-l)/2 



(A31) 



-p/2 



"(9-l)/2 



V p < V < V a 



Va < V 



Eqs.A5-A31 give a complete description of the self- 
absorbed RS spectrum at any time t after shock cross- 
ing. Employing this formulation, Fig.Al shows the time- 
dependent, broadband afterglow spectra of a burst at z = 1 
with our fiducial parameters discussed in §2, where the RS 
contribution has been displayed separately along with the 
total FS plus RS emission. 
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log v [Hz] 

Figure Al. Broadband spectra of a GRB afterglow at 2 = 1 
at selected times after the burst, for the fiducial parameters dis- 
cussed in §2. The solid curves are the total forward shock plus 
reverse shock spectra, while the dashed curves show the reverse 
shock component only for t =10 min, 2 hr and 1 day. 



